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Appl. No. 10/773,371 

Reply to Office Action of September 1, 2005 

REMARKS /ARGUMENT S 

The Examiner is maintaining the rejection of the claims 
primarily over Kobayashi et al. alone or in combination with 
other references. There is also a rejection of claims 17-22 and 
27 over Ichinose; and a double patenting rejection over 
US 6,582,802 to Ohbayashi et al. of claims 17 and 20. 

The Examiner's rejections are based on the Examiner's 
interpretation of the meaning of the claims in view of the 
Examiner's misunderstanding of the term "constant drying rate 
period-" "falling drying rate period." Upon review of the 
Examiner's comments applicants realize that the figure showing 
. the rate of total moisture content and the temperature of the 
target material was misunderstood. This is f in part, because the 
temperature of the target material, as shown in the figure, rises 
during the "falling drying rate period," as explained in the 
specification (see discussion below) . In any case, it is 
requested that the Examiner reconsider all of the rejections, 
including the double-patenting rejection, in view of the 
following further explanation concerning the "constant" and 
"falling" drying rate period as explained below. 
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I. Explanation of the Meaning of "Constant" and "Falling" Drying 
Rate Periods for the Present Invention 

Annexed hereto are two figures to help illustrate the 
intended meaning of the terms in question. Fig, A is the same 
figure previously presented. Fig. B is a new figure showing the 
rate at which moisture is vaporized. 

"A constant drying rate period" is one wherein the same 
amount of moisture which vaporizes per unit of time in the period 
is constant. Therefore, the rate at which the total moisture 
content falls is shown as a constant rate in Fig. A (II). As far 
as a temperature of target material is concerned, the temperature 
of the material remains essentially constant because the heat of 
vaporization of the solvent in the period (II) is the controlling 
factor* In other words, since the solvent is essentially free, 
the temperature of the material is controlled by the evaporation 
of the solvent (heat of vaporization to evaporate the solvent 
keeps the temperatures from rising) . 

"A falling drying rate period" is one in which the amount of 
vaporized moisture reduces in the period. Since the free solvent 
has been driven off in the earlier period, the total moisture 
content is relatively small and tied up in the material. 
Therefore, as a practical matter, for the drying of the present 
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invention, the amount of total moisture content does not greatly 
change as shown in the Fig. A (III). Without the vaporization of 
the solvent, the temperature of the material goes up because the 
material is heated and heat of vaporization does not absorb the 
heat, in the falling drying rate period (III) . 

II. Support in the Specification for the Meaning of "Constant" 
and "Falling" Drvlna Rate : 

Even if the art could interpret the terms differently, the 
specification defines the terms. Attention is directed 
especially to the disclosure in the specification at page 21 et 
seq. which discusses the drying process. Attention is directed 
especially to the discussion in the last paragraph on page 21 
through the second paragraph on page 22. As described on page 21 
starting in line 10, the constant drying rate period occurs when 
the ink-jet recording sheet comprises a relatively large amount 
of water and solvents during the initial drying period. The 
surface temperature of the side having the porous layer is kept 
almost constant due to the fact that relatively free water and 
solvents evaporate while deriving those of latent heat 
vaporization. The constant temperature period is called the 
constant drying rate. 
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The discussion on page 21 continues to discuss the period 
when water and solvents , which interact with the hydrophilic 
binders, are evaporated. It is explained that, since there is no 
free solvent evaporation, the temperature increases (to provide 
energy to release any interaction between the solvent and the 
hydrophilic binders) . As explained in the first few lines on 
page 22, the evaporation rate of water from the surface is 
greater than the movement of water in the coating layer. It is 
during this falling drying rate period in which void formation 
starts. 

In view of the above and the further discussions concerning 
methods of confirming the drying rates, etc. which follow in the 
specification, applicants submit that the meaning of the 
"constant drying rate period" and the "falling drying rate 
period" as previously discussed and as explained above, is fully 
supported by the specification as originally filed. 

III. Art Publications Support A pplicants ' Position 

The Examiner appears to consider that applicants position 
with respect to the meaning of these terms is contrary to the 
meaning known to one of ordinary skill in the art (bottom of 
page 8 and beginning of page 9 of the Office Action) . Therefore, 
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the Examiner considers that applicants 1 arguments are at odds 

with the conventional meanings of the phrases. The Examiner did 

not provide any publications for consideration by applicants 

showing the Examiner's position. However, annexed hereto is a 

publication "Modelling the Drying of Solvent Coatings on 

Continuous Webs," which is consistent with applicants' intended 

meaning. In particular, referring to the summary on page 1, in 

the second paragraph, the second sentence states 

"When the diffusion of solvent to the surface 
is found to be less than the drying rate in 
the constant rate period, then constant rate 
drying is over and falling rate drying 
begins. " 

Also, on page 3 in the second paragraph, there is a similar 
explanation as to when the falling rate begins as noted above in 
the discussion of the present specification disclosure. These 
definitions are essentially the same as the ones applicants 
consider appropriate and disclosed in applicants' specification. 

In view of the above, it is requested that the previous 
arguments distinguishing over the art be reconsidered as being 
fully supported by the present specification and the usual 
knowledge in the art, and that the rejections be withdrawn and a 
favorable action on the merits be provided. Also, 
reconsideration of the double-patenting rejection is requested. 
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If the Examiner is maintaining the rejections on the basis 
that applicants are using terminology which is inconsistent with 
the usual terminology in the art, it is requested that some 
evidence be provided of this usual terminology and also it is 
requested that an explanation as to why applicants' specification 
is not sufficient to support applicants' opsit,ion f be provided. 
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Fig- B 
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(I)a preheating period 

( H > a oonstant drying rate period 

(HI) a falling drying rate period 
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MODELLING THE DRYING OF 
SOLVENT COATINGS ON CONTINUOUS WEBS 



Edgar B. Gutoff 

Consulting Chemical Engineer, 194 Clark Road, Brookline, MA 02146 

617-734-7081 



Prepared for presentation at the 46th Annual Meeting of the Society for Imaging Science 
and Technology, Cambridge, MA, May 9-14, 1993. 
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SUMMARY 

The drying of solvent coatings on continuous webs has been modeled using a 
spreadsheet in Lotus 1-2-3 v. 3. The dryer may have up to seven drying zones where the 
air temperatures, air velocities, and solvent concentrations in the air can be varied 
independently. The web is assumed to be impervious, with evaporation of the solvent from 
one side only, but the drying can be cither single-sided or double-sided. The solvent can 
consist of up to three components. The data base in the spreadsheet now contains 14 
solvents, and more can be easily added. The Wilson equation is used to account for non- 
ideality of the solvent mixtures. Wilson constants for 44 solvent pairs are presently in the 
data base. 

The mode] calculates the rates of evaporation in the constant rate period and in the - 
fall ing rate period. When the diffusion of solvent to the surface is found to be less than the 
drying rate in the constant rate period, then constant rate drying is over and falling rate 
drying begins. 

The heat transfer coefficient at air velocity or plenum pressure must be known or 
estimated. It is then calculated for the particular zone conditions. The web temperatures 
during the constant rate period are found using the equations of Gary and Gutoff . A factor 
used to correct for the reduction of the vapor pressure due to dissolved species can be 
varied to give the correct web temperature in a test run. 

The solvent dlffusivlty In the coating is a function of temperature and of solvent 
concentration, and a three-constant equation based on free volume theory is used to 
estimate it One constant is the activation energy, and that is taken as a value typical for 
diffusion, about 15,000 J/mol (3582 cal/mol). Changing it slightly seems to have little effect 
The other two constants are varied to make the model match a test run. These two 
constants do not seem to be independent and actually only one needs to be varied. 

The coating is divided into a number of slices, usually ten, in which the solvent 
concentrations and the dif fusivities are calculated. The drying rate in the falling rate period 
is equal to the rate of diffusion to the surface. 

When drying air is recirculated the solvent level in the drying air can be calculated 
by a material balance. Otherwise the program assumes that the incoming air has no solvent. 
The concentration of solvent at the surface is assumed to be zero, as if the air has no 
solvent If the air contains solvent then an equilibrium relationship ougjit to be introduced, 
but this relationship is often unknown and is therefore often neglected. This relationship 
could be similar to a published one giving the concentration of water in gelatin in 
equilibrium with air as a function of relative humidity and temperature. 

The spreadsheet consists of four worksheets. Input data is entered on sheet A> the 
calculations are carried out on sheet B, sheet C contains the solvent data base and calculates 
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the average solvent properties, and the macros arc en sheet D. To start the calculations one 
hits Alternated. To change the solvent one hits Alternate-^ A copy of most of sheet A 
is printed out with the temperatures and solvent levels at the end of each zone and also the 
location of the end of the constant rate period. Two graphs, one linear and one semi- 
logarithmic showing the temperatures and solvent levels throughout the dryer, are also 
printed out. 

THE MODEL 

Modelling of solvent drying is divided into three steps: 

- the constant rate period 

- the transition between the constant and falling rate periods 

- the falling rate period, and 

In addition, the heat transfer coefficient at a reference condition must be known. 
All of the modelling, while in principle relatively straight forward, involves a number of 
assumptions which should be explicitly stated. And the modelling for both drying periods 
must be tested to validate the model and to fit constants and correction factors to the 
model. The constants that need to be found are: 

- the heat transfer coefficient at a reference air velocity or plenum pressure. 

- the constants in the diffusivity equation. 

- the vapor pressure factor, . 

With aqueous coating? almost all the drying takes place in the constant rate period, 
and a model just describing the constant rate drying is adequate for many purposes. Such 
a model was developed by Caiy and Gutoff (1991). With solvent coatings most of the 
drying occurs in the falling rate period; in some cases the constant rate period os over 
before the coating enters the dryer. 

In the constant rate period the wet coating behaves as if it were a pool of solvent. 
While the tJ poor remains a/'pool" the evaporation rate will remain constant. The rate of 
evaporation from the "poor of solvent will depend only on the temperature of the solvent 
and on the conditions of the drying air, the rate of flow of the solvent to the surface of the 
coating plays no role at all- In single-sided diying the temperature of the coating will regain 
constant at the wet bulb temperature of the air with respect to the particular solvent. In 
double-sided drying the temperature will be higher. The rate of drying can be found from 
the rate of heat transfer from the air to the coating divided by the heat of vaporization of 
the solvent Modelling the drying mainly involves determining the web temperature. With 
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thin coatings on relatively thin webs the temperature wiH be quite uniform (within a degree 
or two) between the surface of the coating and the bottoxn of the web, except for a fraction 
of a second after conditions, such as air temperatures, change. This assumption of 
temperature uniformity gives us the lumped-parameter model 

At the end of the constant rate period dry patches will appear on the surface of the 
coating and the rate of heat removal by evaporation will no longer equal the rate of heat 
transfer from the air to the coated web. The diffusion of solvent by diffusion to the surface 
by is no longer adequate to keep the surface of the coating wet, the rate of drying 
decreases, and the falling rate period begins. The web temperature increases, rising to the 
temperature of the drying air. When all of the surface appear* dry the rate of drying is 
controlled only by the rate of diffusion to the surface. 

[With porous webs, such as un coated papers, and with porous coatings, such as clay 
coatings with low binder levels, mechanisms other than diffusion will control the evaporation 
rate during falling rate drying. This model only concerns diffusion.] 

BEAT TRANSFER COEFFICIENT 

Before any modelling is done the heat transfer coefficient at a reference plenum 
pressure or reference velocity must be known or measured. It is difficult to calculate the 
heat transfer coefficient, although Martin 1 does give a method for doing this. However, one 
can calculate the heat transfer coefficient from the temperature profile in the dryer of an 
uncoated dry web* This is based on an unsteady state heat balance on the web. If a section 
of the uncoated web is heated by hot air from both sides, then the rate of accumulation of 
energy in the web is equal to the rate energy flows into the web. No energy flows out of 
the web. Thus, following a unit area of web in the dryer gives 

Jr v .o (Tt-Tj Kpe p Jo 

Ww - -H* ln (^kf) (1> 



where « heat transfer coefficient, back side W/m 2 -K 

« heat transfer coefficient, coating side W/m 2 -K 

c p = heat capacity of the web J/kg-K 

= thickness of the web ni 

T t =. air temperature °C or K 

r Wf0 = initial web temperature °C or K 

T w « web temperature °C or K 
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p = web density 
6 = time 

It is better to use a thicker web, perhaps 180 iiro (7 mil), rather than a thin web for 
this test, as it will have a higher heat capacity (h^pc^y 

As reported by Gary and Gutoff 2 the heat transfer coefficient generally varies as the 
0.39 power of the air plenum pressure or the 0.78 power of the air velocity* Thus, if the air 
nozzles are at the same configuration in all the drying zones then the heat transfer 
coefficient is equal to 

h - b W 

or < 2 > 

where h 0 is the heat transfer coefficient at the reference pressure P 0 or the reference 
velocity V 0 . The reference pressure could be several hundred Pascal (about an inchof 
water pressure) and the reference velocity could be about 1000 ro/min (3000 ft/tain). The 
exponent 0.39 could vary by ±0.01 and the exponent 0.78 could vary by ±0.02 without 
significantly affecting the results. In the spreadsheet the reference air velocity is 4000 
ft/min, but can be easily changed. 



kg/m 3 
s 



MODELLING THE CONSTANT RATE 

A good model for constant rate drying of aqueous systems has been published 
recently 3 . The development starts by equating the heat transfer to the web on both sides 
to the mass transfer rate from the surface times the latent heat of evaporation. The 
Chilton-Colburn analogy 3 relates the mass transfer coefficient to the heat transfer 
coefficient. The end result is an implicit equation for the web temperature. 



where P 4 « vapor pressure of the solvent in the drying air on the coating side Pa 

P„ . = vapor pressure of the solvent at the web temperature Pa 

T^toe = temperature of the air on the base side K 

7^ clfi = temperature of the air on the coating side K 
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and the grouping K is 



K « 



c„ +(o „-c M >(^)(£)]^ B (*)(-|) 



where » beat capacity of the air J/kg »K 

Cp, - heat capacity of solvent in the vapor phase J/kg • K 

Af 4 molecular weight of the air dal'tons 

M % «■ molecular weight of the solvent daltons 
Nu* = Lewis number for the solvent in air, (a/Z>) 

a « thermal diffusivity, k/pCp nWs 

Z> » diff usivity of solvent vapors in air m 2 /s 

A = heat of vaporization of the solvent J/kg 

w = atmospheric pressure Pa 

These equations can be solved efficiently for the web temperature by the iterative 
Newton-Raphson procedure. Letting i+J represent the (i+l)th iteration, then 

' .. <« 

where f(T^ is the value of the expression in the Uh iteration, and f f (T^0 h Its derivative 
with respect to T w . Thus 



Here JV is the derivative with respect to the temperature of the vapor pressure of the 
solvent, evaluated at the web temperature. The small variation in the group K with 
temperature is neglected. 

When only single-sided drying is used, these equations will calculate the psychrometric 
chart for the particular solvent, and are similar to the equation derived by Keey 4 and to the 
equation presented by Pakowskf and Mujumdai 5 . 

This model can be used for any solvent or solvent mixture. The vapor pressures of 
pure solvents are found from the Antoine equation 
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log P° (in torr) c A - ^ * - (7) 

The average properties of a mixed solvent can be found from its composition in the 
liquid and in the vapor. For the vapor composition the activity coefficient of each 
component is needed The vapor pressure of one component in a muiticomponent mixture 
is 

where P { = vapor oressure of the hh component Pa 
x\ « mole fraction of the ith component in the liquid 

Yi - activity coefficient of the fth component 

The activity coefficients can be found from one of several equations whose 
parameters depend on the pairwise interactions. Specifically, we use the Wilson equation, 
for which many of the pairwise parameters are given in the massive compilation of 
GmehKng et al.* and are also in the on-line Dortmund data base. Hie equations for the 
activity coefficients using the Wilson equations is given in Appendix A. 

These equations assume that the coating is pure (mixed) solvent. But the coating also 
contains dissolved materials which will lower the vapor pressure. In aqueous coatings we 
have found that this vapor pressure reduction can be neglected with negligible error; in 
solvent coating, where many of the solids are dissolved rather than dispersed, the vapor 
pressure reduction is appreciable. The program includes a vapor pressure factor which can 
be adjusted between zero and one, and which is used to multiply the calculated vapor 
pressures in order to give the correct web temperature in the falling rate period. This may 
be regarded as a fudge factor, for it also will correct for other errors in the constant rate 
calculations. 

The mole fraction in the vapor is found from 

All the physical properties, particularly the latent heat of evaporation, are functions 
of temperature. However, for some solvents this data is not available and, in view of the 
other approximations, we wffl assume that the physical properties are equal to their 
handbook values at about room temperature. 

It has been observed that the web temperature in the constant rate period does not 
change much with distance down the dryer. This implies that one can calculate the initial 
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web temperature, and then assume that the liquid and vapor composition does not change 
further but remains constant at its initial value. This is obviously not so for a mixed solvent, 
for the less volatile species will become concentrated in the coating as drying progress: But 
this assumption greatly simplifies the model and still seems to give reasonably good 
predictions. One could assume that the solvent in the coating is at all times well-mixed and 
then calculate the change in concentration as in a simple distillation, but this obviously not 
be correct One could also assume that evaporation takes place at the surface of the coating 
with no mixing at all of the solvent, but this is also not so, for the surface will be depleted 
with the more volatile component, which- win then diffuse faster to the surface. The 
assumption of a single solvent that does not change composition greatly simplifies the model 
and allows mixed solvents to be handled in the same model as one-component solvents. 



FALLING RATE MODELLING 1/ 
Diffusion 

In the falling rate period the solvent migrates by diffusion to the surface of the 
coating* where its concentration is in equilibrium with the drying air. If the air is solvent- 
free, then the surface concentration is equal to zero. The rate of evaporation is then the 
rate of diffusion to the surface. Hdwcver, the diffusivity of the solvent in the coated layer 
is a function of the solvent concentration, as weH as the temperature. We use a simple two- 
constant model for the concentration effect, based on free volume theory as earlier used by 
Beels and Claes 7 , Including the activation for the temperature effect, the resulting equation 
has three constants 



D m a 



hi] 



(10) 



where A is the pre-exponential factor cm*/s 
B accounts for the concentration effect 
D is the diffusivity! cm 2 /s 

E is the activation energy for diffusion J/mol 
R is the gas constant 8314 J/mol 

/e m1v is the thickness contributed by the solvent in the coating nm 
is the thickness contributed by the solids in the coating Mra 
T is the absolute temperature, K 

The activation energy for diffusion can be assumed to be in the usual range (12,000 
- 21000 J/mol, or 3000 - 5000 cal/mol), say 15,000 J/mol. The constants A and B would 
be found by trying different values in the finished model until the model fits the actual data. 
The final solvent level is decreased when the diffusivity is higher; thus higher values for A 
and lower values of B, which lead to higher diffusivities, predict lower solvent levels. It 
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appears that the constants A and £ affect little else, and that as long as A and B are chosen 
to give the correct final solvent level it does not matter what the individual values are. This 
means that one can choose one constant arbitrarily, and then vary the other constant until 
in a test run the correct final solvent level is predicted. 

Surface Concentration 

The surface concentration of solvent in the coating would be in equilibrium with the 
solvent in the air. If air is not recycled then the concentration of solvent in the air would 
normally be zero, and therefore the surface concentration of solvent in the coating would 
also be zero. This is what we have used in the modeL 

When air is recycled then the solvent level in the drying air can be found from a 
simple material balance. It is easy to add this to the spreadsheet. 

If air is recycled then the surface concentration of solvent would no longer be zero 
and should be determined as a function of the solvent level in the air. The equilibrium 
levels may follow a relationship similar to that followed by moisture in gelatin, as found by 
Gehrmann and Kast 8 . Their results can be recast as 

jr 0,155-0,00075 t zh>0.2 

(in-^) (ID 

X - <0,26 - 0.0015t> <xh>°" rh<0.2 

where rh is the relative humidity 
t is the temperature, °C 

X is the ratio of solvent to solid (here water to dry gelatin) 

One should be aware, however, that some solids are inert, and the ratio X should not 
include insoluble and non-swdlable solids. The determination of the constants for this 
equation is not a simple matter. 

The present model just assumes that the surface concentration of solvent is zero. 
Finite Difference Equations 

To find the solvent concentration of solvent in the coating the coating is divided into 
a number of slices, normally ten, and the diffusion equations are written in finite difference 
form. The finite difference analysis follows standard methods (see, for example, Mickley 
et al except that the solvent Oiff usivities in, and the thicknesses of, each slice differ from 
each other and change with time. The solids in each slice remains constant and equal to the 
total solids coverage divided by the number of slices* 
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Figure 1 illustrates some of the notation for the slices. Position G ia at the base, 
position / at the top of the first slice, position s at the surface. Position m is at the top of 
the mth slice. h m is the thickness of the mth slice, and the diffusivity of solvent in that slice 
is D m . Let us make a material balance on the slice centered at position m; that is, from 
half-way between positions m-2 and m, to halfway between position m and m+1* The rate 
of accumulation is equal to the rate of diffusion into the slice minus the rate of diffusion 
out of the slice, or, considering the time interval from time n to time n+2, 

where C mjn is the concentration of solvent in the mth slice at time n g/cm 3 
AB is the time interval between n and n+Jf s 



Solving for the concentration at the next time interval yields 
n - r , 2A 9 \ C **~* D * - C (3l * £b±) t C *>*.A* ] (13) 

c M c„.»+ 5^VT[ *" n K h *«) h ^ J 

where jD m is evaluated at the mean concentration in the slice, This equation 

is good for all the intermediate slices. For tho bottom slice, the material balance is 

( Mt- )t - 0 - -i^^ir*) (14) 

Therefore 

* 

Co.^ - C 0wD - <C„„ - C lfB ) (15) 

The flow to the surface by diffusion, and therefore the evaporation rate, is just 
Evaporation rate - ° 0 '^ C * j f 16 ^ 

where C g is the surface concentration, and is in equilibrium with the solvent level in the air. 
With dry air, it is zero. 

The. thickness of a slice changes with time because solvent leaves the slice by 
diffusion to the surface. If at time zero we form a number of identical slices, each with s 
g/m* of solids at an average density of PioM , g/cm 3 , and we have W^^m 2 of solvent in 
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the slice, and if there is no change in volume on dissolution of the solids, then the thickness 
of the slice is 



cm 



Ztl* 



CT7 



-2L <» v.my 



(17) 



m 



But the coverage, W^^, is related to the concentration by 



TCt 

», can 3 



(18) 



Replacing in equation 17 by its value from equation 18 gives 

h„ £_ + (i») 



or 



P Bolide 



Pfolv) 



(20) 



There are certain stability requirements for A0 and f^. la equations 5 and 7 the 
concentration at the present time must not have a negative effect on the concentration at 
the future time. This is equivalent to saying that, omitting subscripts, 

-2£f£ < i (21) 



or 



2 2 C22) 



This grouping in equation 22 is called the modulus, and we see that it must be at least 
two for "the computation to be stable. Since the diffusivity is a strong function of 
concentration, we might expect that the modulus would be lowest at the base where the 
concentration of solvent and the diffusivity are highest. On the other hand, because the 
concentration of solvent is highest there the slice thickness would be greatest, which would 
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tend to give a high modulus, It turns out that we have to check the modulus for all the 
slices, and set the lowest modulus equal to two in determining the maximum time interval. 
We may use shorter time intervals. For example, we always want at least ten time intervals 
in each dryer zone. The maximum time interval for stability is recalculated after each time 
period. 

The time interval is directly proportional to the square of the thickness of the slices. 
If the time interval is too short, an excessive amount of computing time may be consumed* 
We start with ten slices, but if the time intervals are too short then the number of slices 
could be decreased. This would involve some changes in the spreadsheet. 



Web Temperature in the Falling Bate Period 

The Web temperature in the falling rate period is found by a heat balance with the 
rate of accumulation of heat in the web plus coating being equal to the input of heat by 
heat transfer from the air (on one or on both sides) minus the heat output by the latent heat 
of the evaporating solvent Thus, in consistent units 

&w ± c p , s )±Z± r iwcr.. w -V'*^ (T -^ fcj y -* J T&f as (23) 

where W i is the coverage of the fth component and the rate of evaporation, AW^AB, is 
taken as positive. This equation allows us to solve for the web temperature after every time 
interval 

END OF THE CONSTANT RATE PERIOD 

At the start, the evaporation rate is calculated for the constant rate period and also 
assuming fallin g rate conditions; that is, assuming a uniform solvent concentration 
everywhere except at the surface. The time at which the calculated falling rate evaporation 
rate becomes less than the constant rate evaporation rate marks the end of the constant rate 
period and the start of the falling rate. 



THE SPREADSHEET 

The spreadsheet consists of four worksheets. Sheet A contains the input data and 
fists some of the important output data. Most of sheet A is printed out, as in Figure 2> 
along with linear and semi-logarithmic graphs of the temperatures and solvent levels in the 
dryer, Figures 3 and 4. Sheet B is where most of the calculations are carried out Sheet C 
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contains the data base for the solvents, and the calculations of the average solvent 
properties arc calculated here- The macros used to run the program are on sheet D. 

All calculations are carried out in SI units., However, the input data are entered on 
sheet A, Figure 2, in whatever units one wishes to use* The data entered. include the wet 
load, the percent solids, the base thickness, the coating speed, and the air velocities and 
temperature. Hie units are converted to SI units when copied to sheet B» Figure 5. If the 
solvent is a mixture the percentages of each component are also entered on sheet A. The 
constants for the diffusivity equation can be changed in sheet A. Sheets A and B can be 
changed so that air plenum pressures are used in place of air velocities. The important 
output data, namely the temperatures and solvent levels at the end of each zone, are 
gathered on sheet A. 

The calculations are dune on sheet B. The far right side, not shown in Figure 5, 
contains a table of all the named variables and their ranges. The average physical 
properties of the solvent are calculated on Sheet C and then copied to the upper left of this 
sheet, see Figure 5, under Constant Rate Properties and Failing Rate Properties. To the right 
of this is the Input Data which is copied from sheet A, and also some Calculated Data. 
Further to the right is a section, shown in Figure 6, where the vapor-liquid equilibrium is 
calculated using the Wilson equation. Antoine constants for the Antoine vapor pressure 
equation and the pair-wise Wilson constants are given here, as copied from Sheet C. For 
a single-component solvent the Wilson constants are not needed. The mole fraction in the 
liquid, x, is found ii) Sheet C and is copied to the right of the Wilson constants. The 
calculations for the mole fraction in the vapor in each zone is carried out just below. The 
Wilson expression Ay (see Appendix A) is given as in the spreadsheet. We assume that 
the composition of a mixed solvent in the coating does not change in the course of drying, 
and therefore the vapor composition, used only to find the vapor pressure, P> and the 
derivative of the vapor pressure with temperature, JP', also does not change In the course 
of drying except for changes in temperature. The justification for this assumption will be 
discussed later. The vapor pressure and its derivative are used to find the web temperature 
in the constant rate regions. 

The constant rate calculations for web temperature are carried out in the table 
beneath Falling Rate Properties, as shown In Figure 5. The Newton-Raphsbn iterative 
procedure is used to solve the implicit equation for the web temperature- The methpd 
converges rapidly, and twelve iterations are normally used. The iterations have converged 
when the numbers in the column headed /(TJ, defined by equation 4, are small* less than 
about Iff 3 . They are usually much less. The heat transfer coefficients arc calculated from 
the air velocities, and the evaporation rate is just the heat transfer rate divided by the latent 
heat of vaporization. 

Underneath this constant rate table is the section headed Location of the End of 
Constant Bate, Figure 7. Here the rate of diffusion to the surface is compared to the 
constant rate evaporation rate. When the diffusion rate is less it means that the constant 
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rate period is over and the calculations would move to a lower section to calculate the 
falling rate drying. The suggested time interval for the constant rate calculations is given 
in the second line below the heading. It must be smaller than the time it takes to traverse 
one-tenth of a zone* because we want ten data points per zone in the summary table. The 
solvent level for each time interval in the constant rate period is calculated in this section. 
When at least one-tenth of a zone is traversed, the data are recorded in a Data Summary 
table. When the distance traversed brings the coating into a new zone, the data for the new 
zone' are brought into this constant rate section, and also the web temperature and the 
solvent level are recorded in sheet A. 

The falling rate calculations are carried out in the table marked Falling Rate 
Calculations, Figure 8, in the section below the Location of End of Constant Rate J Here 
the diffusion calculations are carried out for. each slice. The maximum time interval for the 
falling rate calculations, recorded on the left side below this table, is the lesser of the 
constant rate time interval or the maximum interval for stable iterations, taken as the 
smallest time in the right hand column of the table and recorded below the table in the right 
column. 

When the calculations are concluded the Data Summary table, Figure 9, located to 
the right of the calculations, is completed/ The distances in the dryer can be converted 

* from meters, and the temperatures from Celsius, if desired. Then linear and semi- 
logarithmic graphs of the temperatures and solvent levels versus distance in the dryer, 
Figures 3 and 4* as weH as a copy of most of sheet A, Figure 2, are printed out. 

There is one column set apart on Figure 9, on the left side, and headed To Copy to 
T-web, C. On the actual spread sheet it is to the right of the constant rate calculations. If 
ever the iterative procedure does not converge, and the constant rate table gets filled with 
error messages, ERR, then one has to reset the procedure. One does this by replacing the 
starting point for the web temperatures in the column headed Trweb, C, with reasonable 

• temperatures, such as 25°C (note that a number 25 is just under this little column). Then 
one hits the calculation key, F9, once. This one calculation removes all the error messages. 
Then one should copy the formulas in the column headed To Copy To T-web, C to the 7*- 
weK C col umn, and start again. This procedure is used by the macros when changing 
solvents, and one can do this procedute in an easier manner by starting the macro \Q to re- 
Specify the solvent components, by hitting alternate-^. 

Sheet D contains all the macros and subroutines, with explanations of what they do. 



RUNNING THE SPREADSHEET 

After entering the data start the calculations by running macro \A by hitting 
Alternate-/!. The program takes about five minutes to run. The running time depends on 
the final solvent levels. The lower the solvent levels the longer it takes to run the program. 
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To change solvents hit Alternate-^. Tiic macro \Q shows a screen with a list of 
solvents and solvent pairs, as in Table I, and asks for the first solvent number. It then aslcs 
for the second and third solvent numbers. If the second solvent number is 2ero then a 
single-component solvent is used and its concentration is entered as 100% on sheet -*4. If 
the third solvent number is zero then a two-component solvent is used and the composition 
entered on sheet A is 50% for each. The correct solvent concentratioii then has to be 
inserted. For a three-component solvent the solvent concentrations are left as they had 
been. When two- or three-component solvents are used the correct concentrations have to 
be entered. 

If a mixed solvent is used, the components have to be chosen among those pairs f Or 
which the Wilson constants in the table. A three-component solvent would require three 
solvent pairs, A-B, A-C> and B~C. Sheet C contains data for 14 solvents and 44 solvent 
pairs, out of a possible 91 solvent pairs (sec Table I). It is £ airJy easy to enter new solvents 
and additional data in Sheet C. 
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IS 



NCivIENCLATLtkE 



A constant in Antoine equation . 

A pre-exponential factor in diffusivity equation vails 

Ay^A^ constants in Wilson equation cal/mol 

B constant in Antoine equation . °C 

B constant in diffusivity equation 

Cj heat capacity J/kg-K 

C constant in Antoine equation *C 

C concentration of solvent g/cm 3 

D diffusivity m z /s 

B activation energy for diffusion J/mol 

h thickness m 

h heat transfer coefficient W/m 2 -K 

k thermal conductivity W/m 2 *(K/m) 

K grouping of constants defined by equation 4 Pa/k 

M molecular weight daltons 

A^Le Lewis number, afD 

P pressure or vapor pressure Pa 

J> vapor pressure of solvent in the diving air Pa 

P„ vapor pressure of solvent at the web temperature Pa 

R gas constant J/mol »K or cal/mol 

rh relative humidity - 

s solids in a slice of the coating e/ m 

t temperature °C 

T temperature K 

V air velocity 

W coverage or amount of material g/m 2 

x mole fraction in the liquid 

X ratio of solvent to dissolved solids 

y mole fraction in the vapor 



Greek 

a . thermal diffusivity, k/pc p m 2 /s 

6 time ' s 

X heat of vaporization J/kg 

rr atmospheric pressure Pa 
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Superscripts 

0 pure single-component solvent 

' differentiation with respect to temperature 



Subscripts 

a air 

base base 

ctg coating 

i component / 

m slice m 

n timen 

o initial or reference condition 

s surface 

solv solvent 

solids solids 

w web 
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Appendix A 

Calculation of Activity Coefficients Using the Wilson Equation 

From Gmehling et al. (1977) 

The Wilson equation for the activity coefficient of component / in a multicomponent 
mixture is 



where 



Vj 1 * is the molar volume of pure f 

R is the gas constant 

A 12 « X 12 - An, a pairwise Wilson constant 

= X u - the second pairwise Wilson constant 



cm 3 /mol 
L987 cal/mol- K 
cal/mol 
cal/mol 



PACE 27/39 ■ RCVD AT 10/27/2005 4:35:29 PM [Eastern Daylight Time] ■ SVR:USPTO-EFXRF-6/27 ■ DNIS:2738300 * CSID:+1 212 319 5101 • DURATION (mm-ss):10-22 



OCT. 11. 2005 4:47PM +1-212-319-5101 customer 01933 



NO. 1992 



P. 28 



Table I 

AVAILABLE SOLVENTS AND SOLVENT PAIRS 
Wilson Constants are entered for the listed solvent pairs. 



No Solvent ' Solvent numbers for solvents paired with listed solvent 



1 Acetone 




2 


4 


5 


6 


7 




11 




13 


14 


2 n-Butanol 


1 






5 


6 






11 


12 


13 


14 


3 Cyclohexanone 




• 
















13 


14 


4 Ethanol 


1 






5 


6 


7 


8 






13 


14 


5 Ethyl acetate 


1 


2 


4 




6 


7 




11 




13 


14 


6 Methanol 


1 


2 


4 


5 


* 


7 


8 


10 11 


12 


13 


14 


7 Methyl ethyl 


1 




4 


5 


6 






11 




13 


14 


ketone 
























8 Methyl i-propyl 






4 




6 






11 








ketone 
























9 Methyl n-propyl 














m 


11 




13 




ketdne 
























10 Methylene 










6 




\ 


■ 








chloride 
























li i-Propanol 


1 


2 




5 


6 


7 


8 9 




12 


13 


14 ■ 


12 Tetrahydrofuran 




2 






6 






11 






14 


13 Toluene 


1 


2 3 


4 


5 


6 


7 


9 


11 








14 Water 


1 


2 3 


4 


5 


6 


7 




11 


12 
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D s ,h s 



Dm-1'hm-1 



s-1 



m+1 



m 



m-1 




/ / / / / 
///// Base 



7" 



Figure 1. Finite difference notation for the slices in the coating for the falling rate period. 
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COTCTAX? AND FALLBfe Ml 
Programed by B» B. Cutoff ' 

See BOTES helw 

%t load, 5/1*2 100.D0 SOLVES? 
1 solids 2S.00 Acetone 

Base id Kylar XeCS 
Base thickness, alls S.oo M 

Coating speed, ft/bin 100.00 total 



tost fcaople 1 

v. 3.01 - February 1993 



A, CB»2/s nqoo 

35,00 B 1.6 

30,00 B, J/ftOl 15000 
100,00 



Zone ' Length 
ft 


' oa 
Length 
ft 


us vELoan 

ft/iin 
ctd side bad aids 


Htl 
ctg 
F 


t^dry 
Basa 
p 


• at end 
r 


Soltent 
at end 
g/M 


1 20.0 


20.0 


3000 


3000 


80.0 


86.0 


75.9 


34.200 


2 20.0 


40*0 


3000 


3000 


90.0 


90.0 


87.2 


16.917 


3 20.0 


60.0 


3000 


300O 


100.0 


100.0 


97.8 


8.056 


4 20.0 


BO.O 


3000 


3000 


120.0 


120.0 


UB.5 


3.593 


5 20.0 


100.0 


3000 


3000 


120.0 


120.0 


115.* 


1.681 


6 20.0 


120.0 


3000 


3000 


120,0 


120.0 


U9.6 


0.823 


7 20.0 


140.0 


3000 


3000 


120.0 


120.0 


U9.8 


0.413 



fiie constant rate period Is over 



3 feet into the drjer. 



Figure 2. Output table showing a portion of worksheet A with the input and calculated data. 
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Figure 3, Linear output graph, showing web temperatures and coating solvent levels 
throughout the length of the dryer. 
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1E+03 



1E-01 



Dryer summary 

Test Example 1 




50 100 
Length in dryer, feet 



Figure 4. Semi-logarithmic output graph, otherwise simliar to Figure 3. 
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DRYIRG: COSSET AHD PALLDIG FATS tet Ixuple 1 
Pngruned by e. B. Gatoff V. 3,01 Penary 1993 



b 6 4OD0fpil,W/»tfi-K 100,00 IIFUT UAH . 

lit load, g/a A 2 



CONSTAT RATE PROPERTIES 



Cp vapor, J/Xg-K 
cp alt, J/kg-K 
Bt; vapor, J/g 
lewis no. 

Vap press redaction 

SOLIDS CUATE) 
Density, g/ci*3 
specific hsat 



1311,64 
1005.00 

2.27 
0.50 



1.16 
0.28 



I solids 
Base id 

Tim**. frki-W — 



100.00 
25.00 

lyiar 



•BASB 

Density, g/cT3 1.01 
Specific heat 0.40 
FALLBfG RATE PROFER'TIBS 
Diffusion constants 
A 11000 H-fl 
B 1.60 
E 15000 J/BOl 
at vapor, J/g 721.4 
Specific bt liq .0.558 
Dens flolv, g/a* 0.796 
Niaber slices 10.00 



Coating speed, i/iin 30,48 



calculated daia 

Solvent ctd, g/V3 75.00 

Solids ctd, g/M 25.00 

Solids thicks, m 21.55 

Solvent thickness^ 94.28 

Basevt, g^2 128.27 

total dryer leitgti 30.48 

Mai the in dryer, a 60.0 

Solids tuidt/sllce, m 3.16 



Constant Rate calculations 



Jonelengtn coated SM back side Mry T-dry BvapRate Cm 
i Mr v«l h Air vd h ctg Base g/i«2-« Length 
ftyWn H/V2-K ft/*in l/*M C C l 



1 


6.10 


3000 


79.9 


3000 


79.9 


26.7 


26.7 


4.43 


6.10 


2 


6.10 


3000 


79.9 


3000 


79.9 


32.2 


32,2 


5.05 


12.19 


3 


6.10 


3000 


79.9 


3000 


79.9 


37.8 


37.8 


5.78 


18.29 




6.10 


3000 


79.9 


3000 


79.9 


48.9 


46.9 


7.35 


24.38 


S 


6.10 


3Q0O 


79.9 


3000 


79.9 


46.9 


48.9 


7.35 


30.48 


6 


6.10 


3OO0 


79.9 


3000 


79,9 


48.9 


48.9 ' 


7.35 


36*58 


7 


6.10 


3000 


79.9 


3000 


79.9 


46.9 


48.9 ' 


.7.35 


42.67 



Zone T-¥sl) T-vab K/1C0Q p-dev prab pHreb f(W) f'(M0 



C K »» kPa kPa/l 



1 


8.3 


281.5 0.30688 


0 


5.445 


0.291 


7H5 


52+02 


281.5 


2 


11.4 


284.5 0.30688 


0 


6.176 


0.324 


o&too 


6E4Q2 


284.5 


1 




2B7.0 0.306B8 


0 


7.036 


0.361 


-72-15 


62+02 


287.0 


♦ 


18.5 


291.7 0.30681 


0 


8.862 


0.440 


08400 


71*02 


291.7 


5 


19.5 


291.7 0. 30688 


0 


8.862 


0.440 


OE+00 


72+02 


291.7 


6 


18.5 


291.7 0.30688 


0 


8.862 


0.440 


QE+00 


7E+02 


291.7 


7 


18.5 


291.7 0.30688 


0 


8.862 


0.440 


OE+00 


7B+02 


291.7 



Figure 5. Portion of sheet B showing the data from sheet A converted to SI units, and the 
results of the constant rate calculations. 
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NO. 1 992 P. 35 



Programed by E. B. Gatoff v. J-M February 1993 



b i mtpjujifc* 100,00 ism DMA . 

Nat load, 9/1*2 



COHSTANT RATE PRQPERHK 



fit solvent 
Cp vapor, J/Xg-K 

air, J/kg-K 
Bt. vapo^ J/g 
Levis no. 

Vap press redaction 

SOLIDS COATED 
Density, g/ci A 3 
Specific heat 



50*32 
1311.64 
1005.00 
K0>fil 
2.27 
0.5O 



t solids 
Base id 

Base thickness, (is 
Coating speed, l/iin 



100.00 
25.00 

lyiar 

127.00 
30.48 



1.16 
0.28 



Density, g/ca'3 1.01 
Specific heat 0.40 
FALLING RATE PROPERTIES 
Diffusion constants 
A 11000 1-9 <* A 2/i 
B 1.60 
E 15000 J/lOl 
fit vapor, J/g 711.4 
specific St Uq 0.55B 
Dens solv, g/ca A 0.796 
slices' 10.00 



CALCULATED DATA 

Solvent ctd, g/i*2 75.00 

Bolide ctd, 25.00 

Solids thickness, pa 21.35 

Solvent thickness r pi 94 .21 

Base vt 4 g/s A 2 128.27 

total dryer length 30.48 

Mai the in dryer, e 60 .0 

Solids thfck/Blice, m 2.16 



Constant Me Calculations 



Zone Length coat 


3D SIDE 


BACK SIDE 


Mry 


T-dry 


B?apRate Cm 


i AirVel 


h 


Airvel 


n 


Ctg 


Base 


g/k"2-s Length 


ft/Bin W/s'2-R 


ft/iin Mpl* 


C 


C 


a 


1 6.10 3000 


79.9 


3000 


79.9 


26.7 


26,7 


4.43 6.10 


2 6.10 3000 


79.9 


3000 


79,9 


32.2 


32.2 


5.05 12.19 


3 6.10 3000 


79.9 


3000 


79.9 


37.8 


3?.B 


5.78 18.29 


4 6.10 3000 


79.9 


3000 


78.9 


46.9 


48.9 


7,35 24.38 


5 6.10 3000 


79.9 


3000 


79.9 


48.9 


48.9 


7.35 30.48 


6 6.10 3000 


79i9 


3000 


79.9 


48.9 


48.9 


• 7.35 36.58 


7 6.10 3000 


79,9 


3000 


79.9 


48.9 


4B.9 


* .7.35 42.67 


Zone T-veb T-vab x/iooo 


p-dttf 


p-veb p'-veb f(Mi) 


f'(Mf) T-vsb 
new 


C I 




' XPa 


kPa 


tta/K 




I 


1 8.3 281.5 0.30668 


0< 


5.445 


0.291 


7H5 


52+92 281.5 


2 11.4 264.5 0.30688 


0 


6.176 


0.324 


0E-KJ0 


6S402 284.5 


1 13.9 2H7.0 0.3O6BB 


0 


7.036 


0.361 


-7B-15 


68+02. 267.0 


4 18.5 291.7 0.30688 


0 


6.862 


0.440 




7E+02 291,7 


S 18.5 291.7 0.30688 


0 


8.862 


0.440 


0E+00 


75+02 291.7 


6 18*5 291.7 0.30688 


0 


6.862 


0.440 


os+oo 


7S+02 291.7 


7 18.* 291.7 0.3068B 


0 


8.-B62 


0.440 


QE+00 


TWil 291,7 



Figure 5. Portion of sheet B showing the data from sheet A converted to SI units, and the 
results of the constant rate calculations. 
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NO. 1992 P. 36 



v-iL lutein* Constants tflison Constants 

Solvent vtl cc/iol i B C Hi Vi3 A-i3 I 

icetOTC 35.00 74.05 7.11714 1210,595 K9.664 - -170,1753 -212.9286 0.2B55 

XeOfl 35.00 *40,73 8*08000 1582.300 239.730 594.1757 - 767.5030 0.5175 

EE 30.00 90.17 7.06355 1261.340 224.969 47^87 -217-5920 - 04971 



100.00 



1.0QO0 



tone Bolv ip, Pa t-il Iri2 Iri3 gum P, Pa 1 V 

1 1 14325 1.0000 0.7456 1.B466 1.2517 2621 0.181 

2 6713 0.6255 1,0000 0,5614 U76B 204* 0.375 

3 6062 0.3501 0.6565 1.0000 1.3065 780 0*141 

5445 1.000 291 

2 1 16584 1,0000 0.7432 1.6385 1.1944 2827 0.458 

2 8007 0.6356 l.OOOD 0.5696 1.1758 2436 0.194 

3 7105 0.3533 0.6637 1.0000 1.3036 913 0.148 

6176 1.000 324 

1 18707 1.0000 0,7412 1.8318 1.1926 31B4 0,463 

2 9258 0.6415 1,0000 0,5765 1.1749 2814 0.400 

3 8095 0,3559 0.6615 1.0000 1.3011 1038 0.148 

7036 1.000 361 

4 1 23139 1-0000 0.7277 1.8200 1.1894 1928 0.443 

2 11964 0.6522. 1.0000 0.5889 M734 3532 0,410 

3 10191 0.3607 0.6575 1.0000 1.2568 1303 0,147 > 
— > M« l- m 440 

5 1 23139 1,0000 0.7377 1.6200 1.1294 3928 0.443 

2 11964 0.6522 1.0000 0.5889 1.1734 3632 0.410 

3 10151 0.3607 0.5575 1.0000 1.2968 1302 0.147 
6m — > BB62 1.000 440 



Figure 6, Portion of sheet B showing the vapor pressure calculations for mixed solvents. 
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NO. 1992 P. 37 



LOCMIW OF BSD OF CCBSMfT SUE 



First 










8,1 


W 


8.3 


47.02 teup, F 


0.50 m delta tiw» 


1.20 8 


, for o.l mm 


0 


1 


2 


3 






0.50 


2.00 


1.50 




0.00 


0.25 


0.51 


0.76 


9.65L 2.54 


0.00 


O.OO 


O.OO 


0.00 




1.0B 


1.08 


1.08 


1.08 




0.40 


0,00 


0.00 


O.OO' 




75,00 


72,78 


71,37 


69.15 




75.00 


72,78 


71.37 


6U5 




7.89 


7.66 


7.51 


7,28 




12,01 


11.79 


11.60 


ii.a 




0.65 


0.65 


0.65 


0.64 




13.61 


13.52 


13.45 


13.35 




4.13 


4.08 


4.04 


3.99 


4.43 la CR rata 



tfeb.tflaperature, C 

Delta t, s, suggest . 0.50 tte 

iteration 

Tiie, s 

Location into dryer, t > 
Solvent concn at surface 
adcluusB of half turf layer, 111 
Solvent in half sarf layer, g/b'2 
Total other solvent, g/i"2 
total solvent, g/a A 2 
Solvent in each layer, g/a*2 
ftuctaesa of other layers, yi 
Solvent ooncm in layers, g/cn A 3 
Difftisivity, ;i A 2/5 
Evaporation rate, g/r2-s 



2cme Previous 0 &is me 1 2 is next tone 

Zone length, i 6.10 6.10 

lrmaent of length, a 0.25 
Iterations in constant rate period 1 

length to end of zona 6,10 0.00 to end pro mm 

Bid of Constant Kate, t into dryer 0.76 in bom 1 1 



.Figure 7. The calculations for the end of the constant rate period, end the solvent level in 
the constant rate period. 
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NO. 1992 P. 38 



T>\r ttiv» m*M» WFT ATT fflW 
|*04UUXt« IV1*M » — * 



SUcb Cone aitf Solvent BiM I«rt ltat del 

i (4 bl H C-i,nH 5l*e 

5/<a*3 p» g/a?2 n'i/6 g/q'3 6 

1 0.03 2.24 0.07 941 <>•" 

2 O.03 2.24 O.0? 947 0.03 0.27 

3 0.03 2.23 0.06 9.14 0.03 0.27 

4 0.03 2.23 O.0S 9,08 0.03 0.27 

5 0.02 2.22 0.05 9.00 0.02 0.27 

6 0.02 2.23 0.04 B.B9 0.02 0.27 

7 0.02 2.20 0.03 8.75 0.01 0.28 
B 0.01 2.1B 0.02 6.» 0.01 0.28 

9 0.01 2.17 0.01 0.« 0.01 8.28 

10 0.00 2 4 6 0.00 8.27 0.00 0.28 



StM 22.07 0.41 • . 

Uo.tf slices M.00 

Solids par atfea, g/n*2 2.50 

. iSfcl-ttH *» tine 0.27 

«ie 8 M.26«Bttl» 04.53 

location into dryer, s 41.80 lone lerflth, i «40 

Bvaporation rate, g/»"2-s 0.02 

Som: ftav 6 j Biis= 7 sett w» • 

totetperature.C m tap 48.79 

Mil, ctg side, C 48.9 Bab ta«p ( r 119.82 

Mir, baae side, C 48.9 

h ctg eider, 7 »- 9 
b base side, ir/r2-K 



Figure 8. The falling rate calculations. 
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NO. 1992 P. 39 



MIA StWflRY Iteration 
Count Appro* length tiu Bvaptat T*wb Solvent Delta t length P-veb 
V Seas ■ ; «/S*2-s C 5/S*2 ; ft ? 



71 


0.0 


0.00 


0.00 


4.43 


8.3 75.000 




0.0 


47.0 


n 


0.1 


0.76 


1.50 


1.22 


8.45 


61.47 


0.31 


2.5 


47.2 




0.2 


1.25 


2.46 


0.98 


13.05 


56.35 


0.33 


4.1 


55.5 




0.3 


1.92 


3.77 


0.77 


17.50 


51,76 


0.33 


6.3 


63.5 




0.4 


2.55 


5.09 


0.67 


20.27 


68.07 


0.33 


8.5 


6B.5 




0.5 


3.00 


6,07 


0.62 


21,60 


45.64 


0.33 


10.1 


70.9 




0.6 


3.75 


7,39 


C.56 


22-77 


42.72 


0.13 


12.3 


73,0 




0.7 


4.43 


B.71 


0.52 


23.50 


40.06 


0.33 


14.5 


74.3 


Tt> Copy 


0.8 


4.94 


9.72 


0.49 


23.87 


35.20 


0.34 


16.2 


75.0 


ft 


0.9 


5.62 


11.06 


0.46 


24.22 


35>B7 


0.34 


16.4 


75.6 


J-vabj c 


1.0 


6.14 


12.01 


0.44 


24.41 


34.20 


0.34 


20.1 


75.9 




1.1 


6.03 


13.44 


0.43 


26.97 


32.06 


0.23 


22.4 


80.5 


•273,2 


1.2 


7.33 


14.43 


0.42 


28.37 


30.51 


0.32 


24.1 


83.1 


-273.2 


1.3 


7.98 


15.72 


0.41 


29.39 


28.48 


0.32 


26.2 


-84.9 


•273.2 


1.4 


B.63 


16.99 


0.40 


29.91 


26.52 


0.32 


28,3 


65.6 


-273.2 


1.5 


9.27- 


18.25 


0.33 


30.20 


24.63 


0.32 


30.4 


86.4 


•273.2 


1.6 


9.91 


19.51 


0.37 


30.3B* 


22.82 


0.31 


32.5 


86.7 


•273.2 


1.7 


10.39 


20.46 


0.36 


30.47 


21.52 


0.32 


34.1 


86.8 


"273 .2 


Ii8 


11.03 


a. 72 


0.35 


30.56 


19.89 


042 


36,2 


87.0 


1.9 


11.60 


22. 9B 


0.34 


30.63 


18.35 


0.32 


38.3 


87.1 


2S 


2.0 


12.32 


24.25 


0.33 


30.68 


16.99 


0.32 


40.4 


87.2. 




2.1 


12.80 


25.21 


0.33 


32.54 


16.02 


0.31 


42.0 


90.6 




2.2 


13.43 


26.43 


0.33 


34.55 


14.82 


0.30 


44.0 


94.2 




2 3* 


14.03 


27.61 


0.33 


35.47 


13.7} 


0.29 


46.0 


95,9 




2.4 
fit 


14.77 


29.07 


0.32 


35.99 


12.49 


0.29 


48.5 


96.8 




2 5 


15.36 


30.24 


0.31 


36.19 


11.59 


0.29 


50.4 


97.1 




2.6 


15.95 


31.40 


0.30 


36.31 


10.77 


0.29 


52.3 


97.4 




2.7 


16.54 


32.57 


0.29 


36.39 


10.01 


0.29 


54.3 


97.5 




2.6 


17.H 


33.74 


0.28 


36.46 


9.31 


0.29 


56.2 


97,6 




2.9 


17.73 


34.91 


0.2T 


36.51 


8.66 


0.39 


58.2 


97.7 




3.0 


18.33 


36.09 


0.26 


36.57 


8.07 


0.30 


60.1 


97. B 




3.1 


IB. 93 


37.26 


0.26 


41.55 


7.51 


0.28 


62.1 


106.8 




3.2 


19.61 


38.60 


0.27 


45.19 


6.BB 


0.26 


64.3 


113.3 




3.3 


20*13 


39.62 


0.26 


46.44 


6.41 


0.25 


66.0 


115.6 




3.4 


20.76 


40.87 


0.25 


47.19 


5.88 


0.25 


68.1 


116.9 




3.5 


21.39 


42.11 


0.23 


47.54 


5.*9 


0.25 


70.2 


117.6 




3.6 


-22.02 


43.35 


0.22 


47.73 


4.95 


0.25 


72.3 


117.9 




3,7 


22.65 


44.59 


0.20 


47.85 


4.55 


• 0.25 


74.3 


118.1 




3.B 


23.29 


45.84 


0.15 


47.94 


4.1B 


0.25 


76.4 


.118.3 




3.9 


23.79 


'46.84 


0.16 


48.00 


3.91 


0.25 


78.1 


116.4 




4.0 


24.44 


48.09 


0.17 


48.07 


. 3.59 


0.25 


80.2 


116.5, 




44 


25.07 


49.35 


0.16 


46.12 


3.31 


0.25 


82.3 


116.6 




4.2 


25.71 


50.61 


0.14 


48.16 


3.05. 


0.25 


64.4 


118.7 




4.3 


26.23 


51.63 


0.14 


48.22 


2.B5 


0.25 


86.0 


118.8 




4.4 


26.87 


52.90 


0.13 


48.27 


2.63 


0.25 


88.2 


118.9 




4.5 


27.52 


54.37 


0.12 


48.31 


2.43 


0.25 


90.3 


119.0 




4.$ 


20.17 


55.44 


9,11 


48.36 


2.24 


0.26 


92.4 


119.0 




4.7 


28.69 


56.47 


0.10 


48.39 


2.10 


0.26 


94,1 


119.1 




U 


29.34 


57.75 


0.09 


48.42 


1.9* 


0.26 


96.3 


119.2 



Figure 9. The top portion of the Data Summary. The data for Figures 3 and 4 arc taken 
from this table. 
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